1 



O 
O 

(N 

On: 

CO 

> 

VO 
O 
On 
O 
(N 
O 
-t— > 

Oh! 
<D 
43 



Recent results using all-point quark propagators 

A. Duncan a ,E. Eichten b ,and J. Yoo a 

a Dept. of Physics and Astronomy, University of Pittsburgh, Pittsburgh, PA 15260 
b Theory Group, Fermilab, PO Box 500, Batavia, IL60510 

Pseudofermion methods for extracting all-point quark propagators are reviewed, with special emphasis on 
techniques for reducing or eliminating autocorrelations induced by low eigenmodes of the quark Dirac operator. 
Recent applications, including high statistics evaluations of hadronic current correlators and the pion form factor, 
are also described. 



1. Introduction 

In many cases the extraction of multipoint 
hadronic correlators of physical interest requires 
enormous statistics, which, given the difficulty of 
unquenched simulations in QCD, makes it crit- 
ical to extract the maximum physical content 
of each available configuration. Recently the 
pseudofermion approach Q to the calculation of 
all-point quark propagators has been substan- 
tially improved by the use of mode-shifting, 
effectively eliminating autocorrelation problems 
and allowing the accurate extraction of all-point 
quark propagators for arbitrary sources and sinks. 
In this talk, some applications of these meth- 
ods will be described. We will describe results 
for (a) fits of current correlators to chiral La- 
grangian forms, (b) studies of the pion form- 
factor on 10 3 x20 unquenched lattices (with 200 
MeV pions), and (c) the problem of disconnected 
parts in isoscalar channels. 

2. Methodology 

We begin from a bosonic pseudofermion field 
4>ma with action (m a lattice site, a the spin-color 
index, Q the Wilson or clover operator): 



iJQ j Q<t> 



For fixed background gauge field A, simulating 
the field <fi produces the correlator (where << 



O >> means the average of O relative to the 
measure e~ s ) : 

<< <Pma<f>nb >>S(>p) = {H~ 2 ) ma ,nb 
« <l>ma($ H)nb »S(<P) = (H^ 1 ) ma.nb 

— {Q~ l5)ma,nb 

Unfortunately, low eigenmodes of H frequently 
result in very long autocorrelations for the low 
momentum parts of hadronic correlators. These 
correlations can be removed (or substantially re- 
duced) by shifting the guilty IR modes into the 
UV: define a shifted hermitian Wilson-Dirac oper- 
ator, with Xi,i — 1, . JV the N lowest eigenvalues, 
and vj corresponding eigenvectors: 
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Take = sign(Aj): this evacuates the IR 
part of the quark operator, greatly reducing low 
momentum autocorrelations (see Fig. 1). Then, 
we reconstruct the correct all-point propagator by 
adding in the contribution from shifted modes: 
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Accurate mode-shifted all-point propagators can 
be obtained even for very light quarks, when con- 
jugate gradient inversion dies. Moreover, the 
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Figure 1. Reduction of autocorrelation from 
mode-shifting 
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Figure 3. Fits of pseudoscalar and axial- vector 
correlators from all-point measurements 
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Figure 2. Comparison of single vs all-point prop- 
agators 



all-point propagator represents a great increase 
in statistics over conventional one-point propaga- 
tors:e.g. on a 10 3 x20 lattice, the all-point is the 
average of 20000 (admittedly highly correlated) 
conventional propagators (see Fig. 2). 

3. Applications 

3.1. Fits of current correlators 

The large statistics provided by all-point prop- 
agators allows a high precision calculation of 
hadronic current correlators, which can then be 
fit to the predictions of chiral perturbation the- 
ory, ultimately allowing us to extract parame- 
ters of the chiral Lagrangian directly from these 



correlators. Fits of this sort from pseudoscalar- 
pseudoscalar and axial vector-axial vector corre- 
lators are shown in Fig. 3 (using 800 unquenched 
configurations generated on large coarse lattices 
with the truncated determinant algorithm (TDA) 
HH). The values for M OT and F v extracted from 
these fits are consistent with those obtained by 
conventional fits using smeared-local correlators, 
but we also obtain an accurate reading of the size 
of higher order terms in the chiral expansion. 
From chiral perturbation theory, we have: 

&psps{q 2 ) = 2 Al A + A 3 + A 4 q 2 + A 5 (q 2 ) 2 

q + a 2 

A axax (q 2 ) = -^j- + A 1 + A 2q 2 + A 3 ( q 2 ) 2 
q z + M£ 

The fits to these formulas are shown in Fig. 3: for 
example, for the axial vector case we find that the 
one-loop chiral contribution A 2 is extracted with 
a statistical accuracy of < 1%: 

F n = 0.187 ± 0.011, Ax = 0.8148 ± 0.0038 

A 2 = 0.0777 ± 0.0003, A 3 = -0.00442 ± 0.00002 

3.2. Three-point Functions: the Pion Form 
Factor 

To extract the pion formfactor, we need the 
following 3-point function: 

JtotMf) = E e l ^-^ sm (^)/ sm M 

wxyz 

<i>{z + x, t 2 )-f 5 ^(x, t 2 )jo m (y> h)ip(w, t o h 5 ip(0, to)> 
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Figure 4. Connected and Disconnected contribu- 
tions to Pion Form Factor 

For large time (t = << ii << t 2 ) (using 
optimized smearing functions / sm ): 

J tlt2 (q) -> (M„ + E(q))e- M ^- E ^F n (q 2 ) 

Of course, one needs this convergence before the 
signal disappears! This requires a large ensem- 
ble. The results for F n (q 2 ) using times (t\,t2) = 
(3,6) and (4,8) are shown in Fig. 4, using mode- 
shifted all-point propagators for 65 configurations 
of unquenched 10 3 x20 lattices at M v ~ 200 MeV 
generated by the TDA method. We expect that a 
useful result for the pion form factor will emerge 
from an ensemble of a few hundred configurations, 
which are presently being generated. 

3.3. Disconnected parts and the r\' 

Calculations of the rf mass from isoscalar cor- 
relators suffer from the need to extract discon- 
nected (annihilation) graphs which require all- 
point propagators. Some examples of connected 
and disconnected contributions to the isoscalar 
correlator computed with mode-shifted all-point 
propagators on 5 configurations of unquenched 6 4 
lattices are shown in Fig. 5. The statistical errors 
of the connected and disconnected contributions 
are comparable using all-point propagators, but 
the fluctuations are large from one gauge config- 
uration to the next, requiring a large ensemble be- 
fore an accurate r( mass can be obtained from the 
subtracted amplitude. However, the difficulties 
of disconnected configurations can be avoided en- 
tirely by extracting the rf mass from scalar isovec- 
tor correlators which are dominated at large time 



Figure 5. Connected/Disconnected allpoint cor- 
relators 
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Figure 6. Scalar and Pseudoscalar Correlators for 
Unquenched 6 4 Lattices 

by a r( — it two-body s-wave state. The isovec- 
tor scalar and pseudoscalar correlators extracted 
from 70 fully unquenched 6 4 lattices (a=0.36 F) 
with two degenerate light sea-quarks (M n ~330 
MeV) are shown in Fig. 6: the corresponding rf 
mass is ~ 735 MeV. 
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